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I. INTRODUCTION 


This paper presents an investigation of a class of projection operators 
that arises in the analysis of the approximation of differential equations by 
spectral methods using Chebyshev decomposition. 

Some similar operators have been studied before by Canuto-Quarteroni [1] 
and Maday-Quarteroni [1], but the existing results are not adequate in many 
applications. In fact they forbid analysis for the error of the approximation 
by spectral methods of fourth-order problems and, in several instances, 
second-order problems (see Canuto-Quarteroni [2]). 

We first present some background tools required for our analysis. They 

2 - l/o 

consist of Sobolev spaces relative to the weight w(x) = (1 - x ) L (this 
weight arises in the relations of orthogonality of Chebyshev polynomials). We 
recall and complete results proved by Grisvard [1], [2] concerning 

interpolation theory between these spaces. 

Then we present an analysis of projection operators from these spaces 
into the set of all polynomials of degree lower than N. 

Finally we give an application of the results herein proved to a simple 
test problem. 

We shall give other applications in Maday [1] and shall in a future work 
extend these results to multidimensional domains. Our aim is to apply such 
results to the analysis of the approximation of Navier-Stokes equations by 
spectral methods (see Maday-Metivet [1], [2]). 

For some different notions about projection operators that arise in 
spectal methods, see Tadmor [1]. 
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II. PRELIMINARIES: SOME FUNCTION SPACES 

Notations and Basic Properties 

Let J be an open interval ]a,b[ of ~R (a<b); we consider a weight 
function p(x), continuous over J, satisfying p(x) >. Pq > 0 for any 
x £ J. 

Let us set: 


2 

(2.1) Lp(J) = {<{>:J->-]R|<j> is measurable and ( <(> , <|> ) ^ < +»} , 

equipped with the inner product (<j),^) p = / <j>(x)^(x)p(x)dx. For any integer 
s > 0 we set: 


where: 

( 2 . 2 ) 


H“(J) - )* £ L P (J)| n*n Si0 < -), 


'♦'p,p ' 


this space being equipped with the inner product: 


(2.3) 


s jk. 

«♦.♦», 0 - i (H - ^ • 

S,P k=0 dx K dx k p 


Clearly, one has the equality: 


L^(I) = H°(I). 
P P 


For any real s 0, noninteger, H S (J) 

S o+l _ 

the space H (J) and H (J), where s 
P P 5 


is defined by interpolation between 
represents the integral part of s. 
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The method of interpolation can be the complex one, the operator's domain one 

or the trace one (see Lions-Magenes [1] for more details). Besides we define 

s s 

H 0 p(J) as being the closure of P(J) in Hp(J). When p = 1 these spaces 

s s 

are the usual Sobolev spaces denoted by H (J) and Hq(J) respectively. For 
the application to spectral methods we are mostly interested in those spaces 
when J = I = ]-!,+![ and p(x) = o)(x) = - • Let us recall some 


A 


2 

- x 


results proved in Grisvard [1], [2] valuable for J = I, p = a) and for 

j = ]o,i[, p = -L . 

/ X 


THEOREM 2.1 (Grisvard [1]): 

i) For any real s > 0, s (t IN + V 4 we have: 


(2.4) 


h o,p ( j) ■ h o!p (j) ] [s 


ii) For any integral n we have: 


(2.4') 




n+ V 4 


iii) For any real s>0, s^lN+ V 2 


H®(J) C C m (J), 

the space of continuous mapping defined over J whose derivative of order 
m are continuous over J, with m = s - y • 

The trace application defined from C (J) into IR? n : 
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u %1>, U(l), £<1>,- f>) 


can be extended to a continuous mapping from ^ +e (I) onto lR? n for any 


e > 0. 


iv) For any real V 4 < s < 5/4, (J) coincide with the subspace of 

u , p 

Hp(J) of functions vanishing at the real boundaries of J. 
v) For any real s > V 2 > H®(J) is an algebra. 


THEOREM 2.2 (Grisvard [2]): 

For any 0£q<s<p,H^(J) 
imbedding: 


satisfies the following double topological 


(2.5) 


[h’W). H p( J >]e,i CH p (J)c 


with 0 = ^ ~ P , and the notation holds for the real interpolation (see Lions 
and Peetre [1]). 

The two following results can be found in Canuto-Quarteroni [2] and 
Maday-Quarteroni [1]. 


THEOREM 2.3: 

i) For any real s > V4 > H S (I) C H® ^ (I). 
ii) For any 0 _< r < s, the imbedding H^(l) ^ H^(I) 


is compact. 


In the next section we shall generalize the results (2.4) and (2.4'). 
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Some New Results About Interpolation Between 

This section is devoted to the proof of the following: 

THEOREM 2.4: For any 0 < q < s < p not in H + V 4 we have: 


[H^ (I), h£ (1)1 = H® (I). 

L 0,a) 0,0) J r s~p -I 0,(i) 

Lq-p - 1 


This theorem is a consequence of the two following lemma: 


LEMMA 2.1: For any integer p < n, we have: 


u £ H S,„ (I) 


4 u £ (X) 

dx p „4(n-pHl 


PROOF: It is an easy matter to check that this result is a consequence of 


(2.6) 




0, 


/x 


dx* (-z) 
/x 


(we shift the difficulties at ±1 onto 0). So let u be in H n .(0,1); 

0> ~P 

from Theorem 2.1 (point iv) we have, for any 0 < p < n: vx 


d^u. 


dx* 


0 ) = 0 , 


x jP+1 

1 


d^Uv 


0 dx* 


hence 
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Besides, from Lemma 6.2.1 of Neifas [1] we have, for any a < 1 and any v 
1 2 a 

such that j v (x)x dx < ® : 

0 


1 x 


(2.7) / (/ |v(x)|) 2 x° 2 dx ^ (yM J |v(x)| 2 x° dx, 

r\ r\ 1 Ot « 


0 0 


1 JF 

taking then o = - y - 2(n - (p+1)) and v we obtain: 

dx p 


1 '^) 2 x“ V 2 ' 2< "" p) dx < C / (i^) 2 x~ 1/2 ‘ 2 <--<>> +1 » dx, 


i 

0 dx P 


0 dx* 


and (2.6) holds by induction over p. 


LEMMA 2.2: For any integer n > 0, the mapping u \ ► uw /2 is an 


homeoraorphism from ^(1) onto Hq(I). 


PROOF : Here again we prove the result for the weight — , say: 

/x 


(2.8) u 


k ^ is an homeoraorphism from (0,1) onto H^(0,1), 

U j 1 u 

/x 


Let <j> € P(0,1), then, for 0 < m < n : 


A**- 1 *). " 

m T *-• m r\ m— t-k 


n m p m-p 

p=0 dx r dx r 


= l C P D P d^ x -l/4~(ni-p) 

“ m m Tk * 


m m - p 
p=0 dx 


with D p 
m 


["Va" (m-p+l)lD p and D m = 1. From Lemma 2.1 we then get: 
J m m 
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'^+ X ' 1/4 >'0,1 i i C, *'n,l//S 1 

dx 


summing up these estimates for 0 m n we derive: 


(2.9) 




Inversely, let us prove that, for any <J> €#(!): 


( 2 . 10 ) 


"^ X 4 "n,l//^< C ^"n,l- 


From Hardy's inequality (Lemma 2.5.1 of Ne£as [1]) we derive by induction 
that, for any 0 < p < m < n: 


(2.11) 


. < C II 4, II . , 

dx? - "- 1 


besides : 


>X ^ ) 


= l c P - m - P (x 1/4 ) 

“ m n m— -n 


n m J P a m ~P 
p=0 dx dx 


= l c P D 'pd!i x VM 

“ m m n 7 


m m , p 
p=0 dx r 


with: 


K V = [ 1/ 4" (® - P + D] 0 ; 1 * 1 and D' m = 1, 


Then using (2.11) we get: 
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d m (4>x 


V4 


dx 


m 


£L (0,1), 

l//x 


and (2.10) is derived by summing up these results for 0 _< ro < n. We can now 
achieve (2.8) as a consequence of (2.9) and (2.10). 


We can now prove the main result of this section. 

PROOF OF THEOREM 2.4: From (2.4) and Lemma 2.2 we deduce that the 

V? 

mapping u | ► urn is an homeomorphism from Hq ^(1) onto H®(I) for any 

s >_ 0 not in I + V/ini + V 2 (see Lions-Magenes [1] for more details about 
the properties of spaces of interpolation). 


Let us recall that, for any q s _< p not in IN + V 2 we have (see Lions- 
Magenes [ 1 ] ) : 

< 2 -‘2> Hq(I) - [H’(I), HP(I)] s _ q . 

p-q 

From the previous homeomorphism we deduce that, for any q < s < p not in 

{ IN + V2 } U { IN + V4 } : 


(2.13) 



p-q 

Let us remark now that the values of p, q, s in IN + V 2 have only been 
excluded due to (2*12), these values can now be recovered thanks to the 
reiteration theorem (Theorem 1.6.1 of Lions-Magenes [1]). 
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III. APPROXIMATION RESULTS OF PROJECTION OPERATOR IN WEIGHTED SOBOLEV SPACES 

r 

The previous theorem leads us to define over H (I) a new scalar 

U 

product. Indeed, for p not in IN + V4UIN + 1/8, (I) can be seen as the 

U , (1) 

interpolate V2 between L^(I) and Hq^O) and for p in IN + 1/8, Hq ^(1) 

2 3n 

can be seen as the interpolate 1/3 between L (I) and H n p (I). 

0) U , 03 

If we consider the domain operator interpolation, this find expression in 
the existence of a self adjoint operator such that: 

3 3 

* if r € IN + 1/8, the domain D(A^) of the operator A f in 

L^(I) is H^ r (I) if r 4 IN + 1/8, the domain of D(A^) of 
to 0 r 

the operator A^ in L^(I) is H^ r (I). 

K r 03 0,03 

* The domain D(A ) of the operator A in L^(I) is (I) if 

r r oi u , oi 

r 4 IN + Va and is included in (I) if r €1 N + Va • 

u , oi 

Moreover: 

(3.1) (u,v) I ► (((u,v))) = (A u, A v) , 

r y a) r r 10 

is a scalar product whose associated norm is equivalent to the one defined in 

(2.2) if r i M + V4 . 


r r 

Let us define now P r ^ as the projection operator from Hq ^(1) over 
with respect to the previous scalar product with: 


S N ' S N r)S 0,J 1) ’ 


S N = { <{> defined over l| <(> is a polynomial of degree _< N}. 
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LEMMA 3.1: Let 0 _< v < r < a with a 4 IN + V 4 we have, for any 

(3.2) « - P r , N *» v>(iJ < CN V "° »*<i a>aj . 

REMARK 3.1: The case v = r = 0 has been studied in Canuto-Quarteroni 

[1], the case 0 _< v £ r = 1 has been looked at in Maday-Quarteroni [1] (note 
that the dependence of the constant is then C(a) = C*(a!)). Moreover it is 
proved that no optimal bound was possible for Hq ^(1) norms with v > r. 
Indeed, for example: 

(3-4) '♦ - p 0 ,N < Cn2 '"° 

It is often necessary (see Canuto-Quarteroni [2], Maday-Metivet [2], Maday 
[1], and (4.14) ) to obtain optimal results in higher norms. 


PROOF: We shall only consider the case r 4 IN + 1/8 for simplicity. 

The proof is divided in two stages 

i) We first prove (3.2) by induction over r in IN . So, let us 
assume that (3.2) is true for s < r in IN ; let € Hq ^(1); then 

€ *C (I) and P r-l,N-l ( +x ) € S N-1 * Moreover if = < f ,(1) = 0 we 

have: 


« * P r-l,N-l<V< t)dt ' [ t [ P r-l,N-l ( V - 


From the Cauchy-Schwarz inequality we derive: 
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l/ 2(( l at-) ^2 


l«l < (/ (Vi )N -l ( ^ " * x > (Oc(t)dt) '2^ (»<t)) 


< cup ,„,(♦)-* i n ; 

— r-l,N-l Y x Y x 0,o)’ 


hence, from the induction hypothesis: 


(3.5) 


I ot I < cn 1- ° n<j. i . 

1 1 — Y x o-l, ( 


Finally we have: 


V*> -/ --f 


ad - tV ’ 1 I., c „r 
dt £ S N* 

/ a - x 2 ) r_1 dx 


Due to the Poincare-like inequality, the polynomial satisfies the following: 


Vr,™ < '<♦ - Vx'r-l,, 


the induction hypothesis, and (3.5) gives us: 


R N l r u < C(N (r N (1 a) ) «<f> x ll a _ 1 ( 


< CN r 0 ld»I 
— Y a,w 


From the equivalence of the norms H II and 

r,co 


, and the identity: 


P M <f> I I I = inf I I I <b - <(> M II I > 
r,N TIM r,aj c r v N 11 1 r ,o)’ 


% C S N 
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we obtain for any (p 


in H°(I)nH„ (I): 

CJ U ,03 


(3.6) 


p . T <(>||| < CN r-CT n<j)n 

r,N T 1 11 r,(o — Y o,b) 


Besides, since the operator is selfadjoint, we have: 


n<f> - p „ <j>n n 

Y r,N y 0,o) 


inf 

* £ I 4 <1) 


( » - p r-N 


lipl 


0,co 


inf 

* £ I 4< 1 > 


( y ♦ - p r ,„ ♦>, a; 1 ♦>, 


nipii 


0,oo 


From (3.1) we then get: 


( (( * ~ P r,N ♦» A r 2 ^^r.to 


(3.7) 


P xt <Mn 
r,N r 0,oo 


inf 


nipn 


0,a) 


By definition of P r N we have, for any ip in L^(I): 


(«♦ - P r,N P r,N^ 2 ♦)»),,„ * »i 


hence 


- P 


r,N 9 0,o3 


* € 


inf 

l »< i > 


«u - p r-H ♦.(a ; 2 ♦) - p r „(a ; 2 ♦)») 


(0 


Hip II 


0,to 


- P 


r,N 


r,to 


* € 


inf 


i i i ti - p r , N (A; 2 


r , o) 


n T^n 


0,03 
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Due to (3. 


Now, from 

we derive 

Due to (2 
(3.8) 


6) we then derive: 


OA 2 i|>ll 


Id) - P „ AH,, < CN r 0 ll<f>ll N r inf -irT„- 2r,o) 
Y r,N Y 0 ,oj — Y o 0) 9 Hi|)H n 

. 1 . /■ u » u 


* £ k,«> 


< CN 0 HAD inf 

— T a a) ^ 

* € L 0) (I) 


< cn 0 ii d> n 

— 0,0) 


,A r( A r 2 0 >o,, 


!^I1 


0,a) 


the two estimates, valuable for any $ € H a (I)OH^ (I): 

0) u , 0) 


p <f.ii < CN r 0 hah 
r,N Y r,o) — Y o 


,-a 


n<t> - p „ Ai n < cn H(j>n 

Y r,N Y 0,o) — Y a,0) 


that for any 0 € ] 0 , 1 [ : 


H - p r N 2 r , < 

[\, (I) - H o,„ CI) ]e 


< CN 0r_O II < 


G,0) 


,4) and (2.4") we deduce that, for any 0 < v < r: 


i<t> - p M 4>n < cn v 0 n 4>n 

Y r,N T v,o) — Y o,o) 


ii) Let us now prove (3.3) for nonintegral values of r. Let 
€ D(I), from step (i) we know that, for any o >_ r + 1 , o (( IN + V 4 
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inf 


*N € S f ' 


- <j> M n < cn 1 ^ 1 a n<>n 

^ r+T a) ~ 


inf 

i € oT+i 

*n € s n 


- <k T n < cn a n<t>n 

N — — r a.U) 

r,u ’ 


(see (3.8) with v = r + 1 and v = r respectively). These two estimates 
are equivalent to the following one: 


n <j>n 


< CH 1 * 1 - 0 I 


C (i)/s r 


0 , 0 ) 




< CN r a hah 

— T 0,0) 


Due to the interpolation of quotient spaces (see Lions-Magenes [1] Lemma 
1.13.2) we have, for any 0 € ] 0 , 1 [ : 


n<f>n 




< CN r+e a n<j>n 


0 , 0 ) 


From (2.4), (2.4') we deduce (we take 8 = r - r): 


«<!>[ 


< CN r ° n< 




I , 
0 , 0 )’ 
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so that, for any 


a > r + 1 a 4 IN + V 4 : 


(3.9) 



d>n < CN r a 
Y r,a) 


n cf) n 


a w 


By definition we know that, for any 


r 4 IN + V 4 • 


n<|> 


- p Ji 

r,N r 


< cn<f>« 

,o) — r, 


a) 


hence, ( 3 . 9 ) holds for any a >_ v, a 4 IN + V 4 and any r £ 1 . The 

estimates in lower order norms are obtained following the same lines as in i). 


We are now interested in the approximation of the spaces H r (I)OH^ (I) for 

(l) , (l) 

s€lN,0_<s<(r by polynomials therein contained. 

For simplicity of expositure we shall consider the case s = 0. From 
point iii) of Theorem 2.4 we can easily exhibit for any $ € H^(I) a 

polynomial (J)q of degree _< 2 r - 1 such that: $ - 4 >q £ Hq ^( 1 ) anc * f° r an Y 

real p: 


(3.10) 


H(b n ii < c ii 4> ii 

Y 0 p , to — Y 


r ,u) # 


Due to the previous lemma, we have, for any 


0<v£r_<a 4 IN + V 4 : 


H* - 4> Q ) - (P r>N (4» " V^'v.a, 1 CN 


v-a 


n - <j> n ii , 

Y ^0 a or 


and (3.10) then implies: 


'♦ - ( *0 + P r,N ( * " V )n v,a,l 


CN 


v-a 


n<j>« 


a a) 
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For N large enough (more precisely N > 2r - 1) we then get the existence of 

an operator P from H r (I) onto S w such that: 

0) w 


I* " P W«v lU < CN 


This estimate provides an answer to our question in the case s = 0. The same 

r s s 

proof can be done to build an operator from H^(I) Hq ^(1) onto 
satisfying analogous bounds. This leads us to state the main theorem of this 
paper: 


THEOREM 3.1: Let (v,r,o) € IR , a i U + V4, s€ IN, 0 _< s £ r, 

s 0 r g 

0 < v < r < a. There exists an operator II ’ from H (I) onto 

— — — r r,N ( 1 ) 

sf. such that, for any <f> e H a (I)P\H^ (I) we have: 

N 0) U ,03 


I < CN V a Hd.ll 

V,01 — T 0,03 


IV. AN APPLICATION 

Definition of the Problem 

In order to explain how the previous results can be applied, we shall 
study an approximation of the very simple problem: 

Find i|. defined over I such that: 
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(4.1) 


djt _ f 

dx 4 ' 




over I, 


at ± 1 . 


(This problem provides a first step for the analysis of Stokes and Navier— 

Stokes problems in the \|>-formulation; see Maday-Metivet [1], [2].) Let us 

-2 • 

define H (I) as follows: 

0 ) 


H w 2 (I) = {f e P'(i)| 3 g e L 2 (I): f=-^}. 

dx 


We now want to prove the following: 


THEOREM 4.1: Let f € (I); then there exists one and only one 

solution ij> to the problem (4.1) in the space H 2 (1). 

0,m 


This theorem is a very simple consequence of Lax Milgram lemma and the two 
following lemmas. 


LEMMA 4.1: There exist two positive constants 6 ^ and 62 such that, 

2 

for any d> in H- (I): 

0,o) 


{ i *1 / £- “ 5 < S 2 I iHf - 

1 I 1 dx 


This lemma is a corollary of Lemma 2.1. 
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LEMMA 4.2: There exist 3 positive constants a, 8 y such that for any 

(♦.+) e H ^ I ) xH o,« u): 

( 4 . 2 ) <♦„.<*•)„> > 

< 4 - 3 > < S'+kx'o,-- 

< 4 - 4 > 


PROOF : 

i) We first note that (4.3) is an easy consequence of the previous 
lemma, and is an equivalent to the Poincare inequality, 
ii) Next, we get the following equalities, for any e V(I): 

[ +xx ( * u) xx dx = { +L “ dx + 2 { *xx *x w x dx + { *xx *“xx dx 

= J 0 )dx + J (^ 2 ) x 0) x dx + J (* x ) x dx 

= J 4< 2 x mdx - / * 2 m xx dx - J * x (+“x*)x dX 

= J <j) xx rndx - 2/ ^ o) xx dx + V 2 J f <» xxxx dx 
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let us note that: 


w xx = (1 + 2 x 2 )(o 5 , 


“xxxx = ( 9 + 72x2 + 2 4x 4 )w 9 ; 


hence: 

(4.5) 


( *xx ( *“ ) xx dx = / < x dx " ^ ^ (1 + 2x2)o)5 dx +1 /2 / 4> 2 (9 + 72x 2 + 24 x 4 )co 9 dx. 


j J- xx J ' X 


Besides, let us set: 


^ - / (^ xx w + ^ x ^ x w + (2x 2 + 10 V 5 ) 2 a) *dx; 


P is ^ 0 and, an easy calculation gives: 


p = / *j ac <*“> xx " 2* 1 0~ 2 //</-/ (5,78x 2 + 0,4839)/ a) 9 , 


so that 


/ 4 » 5 < 5 °/ ♦«<♦«>«. 


from which we derive: 


/ * 2 d + 2 * 2 )/ < 150 / * xx (*o>) xx . 


Using that inequality in (4.5) we obtain a > 0 such that: 
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{ <, XX ( *“ ) XX > “ 6 { * xx “• 


and (4.2) is an consequence of (4.3). 

iii) Finally, let us note that for any (<}>,40 in H 2 (I)xH 2 (I) we have: 

(a) U , (0 


(4.6) " <* xx .1> xx >» + 2J ♦„ t x “ x + i * xx ♦». 


XX 


The following inequality is simple: 


(4.7) 


I (<f> ,tj> ) I < n<i> I. nip il . 

1 Y xx’ r xx 0 ) 1 — Y xx 0 , 0 ) Y xx 0 , 0 ) 


Let us examine the second terra: 


1/ 'f'xx +X “xl = ^ ♦xx^x “x “ ^“1 < 1/ +xx “I 1/2 ^ *x “x “ *1 1/2 ; 


J XX ■ ’X X 


. 2-125 J . 

since to to = x to , we derive from Lemma 4.1 that: 

X 


(4.8) 


1 1 ^xx ^x °x I — C I ^xx 1 0 , 0 ) I ^xx 1 0 , 


o> 


We obtain, in a similar way: 


(4.9) 


| / <f> *o) | < cU L U L , 

i xx xx 1 — 1 XX ’0,0) 1 XX ’0,0) 


so that (4.4) is a consequence of (4.6)— (4.9) . 



- 21 - 


PROOF OF THEOREM 4.1: Let f = g vv be in fT 2 (I). Problem (4.1) is 

CO 

equivalent to the following: 


(4.10) 


Find in V = ^(1) such that, for any <j) in V: 


{ ♦xx^xx ■ { S<« xx - 


The bilinear form a defined by: For any (x,<|>) in V^: 


a(<f>> X ) =* / ♦xx^^xx’ 

is continuous and elliptic over V (see Lemma 4.2), and Lax-Milgram lemma 
gives the existence and uniqueness of a solution of (4.10) hence of (4.1). 


Approximation of Problem 4.1 

We are interested in approximating the solution of (4.1) by a polynomial 
of degree N. We use a Galerkin method approach known as Spectral Method 
(see Gottlieb-Orszag [1] for more details); hence from (4.10) we derive an 
approximate problem: 


2 

Find i|/ N in = S N such that, for any 4> in 


{ W* u) xx ’ { g( *" ) xx- 


(4.12) 
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From Lemma 4.2 we know that problem (4.12) is wellposed in the sense that 
there exists one and only one solution. Moreover, we derive from (4.10) and 

(4.12) : 

a(if> - iJ> N ,<{>) = 0 for any <J> in V.., 

2 0 

so that (remind H ’ T \!> € V, T ): 

2 ,N T N 

(4.13) a(i|> - - t|> n ) = a(f - ^ . 

Due to Theorem 3.1 and Lemma 4.2, we then obtain the following 


THEOREM 4.2: There exists one and only one solution to problem 


(4.12); moreover it verifies, as soon as e H a (I)(~)H_ (I): 

(l) u , CO 


(4.14) 


♦ - 9 < CN 

T T N Z,(o — 


2-c 


a, w 


REMARK 4.1 : The previous estimate is an optimal one in the sense that no 

2 

polynomial of is asymptotically nearer from the solution than the 

solution of the approximate problem. 


REMARK 4.2 : The previous theorem will be extended in a future paper 
where we shall consider a pseudospectral method (much more efficient from a 
computational point of view) for approximating a one-dimensional fourth order 


equation. 
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